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A low-momentum nucleon-nucleon (NN) potential Vi ow —k is derived from modern 
realistic NN potentials by integrating out the high-momentum modes. The smooth 
Vlow—k may be used as input for nuclear structure calculations instead of the usual 
Brueckner G matrix. Such an approach eliminates the energy dependence one 
finds in the G-matrix approach, allowing this interaction to be used directly in 
Hartree-Fock calculations. Bulk properties of ie O have been calculated starting 
from different NN potentials. Our results, obtained including up to second order 
contributions in the Goldstone expansion, are in good agreement with experiment. 



1. Introduction 

A central problem in nuclear theory has long been the calculation of nu- 
clear structure properties starting from realistic nucleon-nucleon (NN) po- 
tentials. There are various models for these potentials, such as the Bonn 1 
and CD-Bonn, 2 Paris, 3 Nijmegen, 4 and the new Idaho potential. 5 They all 
describe the observed deuteron and NN scattering data very accurately. 
However, owing to their strong repulsive core, none of them can be used 
directly in nuclear structure. To overcome this difficulty, the Brueckner 
G matrix has traditionally been the starting point, but, as is well known, 
its energy dependence is an undesirable feature, in particular when dealing 
with Hartree-Fock calculations. 

In this paper, we make use of a different approach to renormalize the 
short-range repulsion of realistic NN potentials, which is motivated by 
the recent applications of effective field theory (EFT) and renormalization 
group (RG) to low-energy nuclear systems. 6 ' 7,8 ' 9, 10 
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A fundamental theme of the RG-EFT approach is that the details of 
the short- distance dynamics are irrelevant for physics in the infrared re- 
gion. One can therefore have infinite theories that differ substantially at 
short distance, but still give the same low-energy physics. 6 ' 10 In nuclear 
physics, various meson models for Vnn, sharing the same one-pion tail, 
give the same phase shifts and deuteron binding energy even if they differ 
significantly in the treatment of the shorter distance pieces. 

Motivated by these considerations, we derive a low-momentum NN 
potential Viow-k 11 '^ 2 by integrating out the high-momentum components 
of Vnn in the sense of the RG, 6 ' 10 and use this smooth potential in nuclear 
structure. 

In this paper, we calculate the bulk properties of 16 in the framework 
of the Hartree-Fock theory, using VJ OM ,_fe's derived from different NN po- 
tentials. The paper is organized as follows. In section 2 we describe our 
method for carrying out the high-momentum integration. Section 3 is de- 
voted to the comparison of our results with experiment. A summary of our 
conclusions is given in the last section. 

2. Outline of Calculations 

The first step in our approach is to integrate out the high-momentum com- 
ponents of Vnn- According to the general definition of a renormalization 
group transformation, the decimation must be such that low-energy ob- 
servables calculated in the full theory are preserved exactly by the effective 
theory. Once the relevant low-energy modes are identified, all remaining 
modes or states have to be integrated out. 

For the nuclcon-nuclcon problem in vacuum, we require that the 
deuteron binding energy, the low-energy phase shifts, and the low- 
momentum half-on-shell T matrix calculated from Vn n must be reproduced 
by Vi ow - k - 

The full-space nuclear Schrodinger equation may be written as 

*„ = H = H + Vnn, (1) 

where Hq is the unperturbed Hamiltonian, in this case the kinetic energy. 
The above equation can be reduced to a model-space one of the form 

PH cS P*^ = £„Ptt„; H eS = H + V cS , (2) 

where P denotes the model-space, which is defined by momentum k < 
k C ut = A, k being the relative momentum and k cu t a cut-off momentum. 
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The half-on-shell T matrix of Vnn is defined as 

T(k',k,k 2 ) = V NN (k',k) + 

, q 2 dqV NN (k', q) fc2 _ ± - ^ + T(q, k, k 2 ) , (3) 

and the corresponding model-space T matrix given by Vi ow -k is 
T low ^ k (p',p,p 2 ) = Vi ow - k {p',p) + 

A 1 

q 2 dqVi ow - k (p', 9 , ., , 

q p z — q z + i0 + 

Note that for Ti ow ^k the intermediate states are integrated up to A. 

It is required that, for p and p' both belonging to P (p,p' < A), 
T(p\p,p 2 ) = Ti ow -k(p\p,P 2 )- In Rcfs. 11,12 it has been shown that the 
above requirements are satisfied when Vi ow -k is given by the folded-diagram 
series 

Vlow — k Q-Q' J Q + Q' JqJ Q-Qi JqJqJ Q + ... , (5) 

where Q is an irreducible vertex function, in the sense that its intermediate 
states must be outside the model space P. The integral sign represents a 
generalized folding operation, 13 and Q' is obtained from Q by removing 
terms of first order in the interaction Vnn- 

The above Vi ow -k can be calculated by means of iterative techniques. 
In Rcf. 14 some iteration methods have been proposed which are suitable 
for non-degenerate model spaces. The method that we use here, which we 
refer to as Andreozzi-Lee-Suzuki (ALS) method, consists in constructing 
solutions of the Lee-Suzuki type. The ALS method converges to the lowest 
(in energy) d states of H, d being the dimension of the P space. We have 
found this method to be very convenient for our calculations. 

Since the Vi ow -k obtained by this technique is non-hermitian, we have 
hcrmitized it by using the simple and numerically convenient procedure 
suggested in Ref. 14. 

We have verified that the deuteron binding energy and the phase shifts 
up to the cut-off momentum A are preserved by Vi ow —k- By way of illustra- 
tion, we compare in Fig. 1 some phase shifts calculated using a cut-off 
momentum A = 2 fm _1 (squares), with those obtained using the full Idaho 
potential 5 (continuous line). It can be seen that the phase shifts from the 
full Vnn are well reproduced by those obtained from V ow -k up to the value 
of Ei a b corresponding to the cut-off momentum. 
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Figure 1. Comparison of 1 5q phase shifts given by Vjvjv an d Vlo 
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This evidences the physical equivalence of V\ ow -k and Vnn in the sense 
of RG. 

An important issue is what value one should use for A. Guided by 
general EFT arguments, the minimum value for A must be large enough 
so that Viow-k explicitly contains the necessary degrees of freedom for the 
physical system under investigation. For nuclear structure calculations A 
turns out to be around 2 fm -1 . This is consistent with the following consid- 
erations. Most AW potentials are constructed to fit empirical phase shifts 
up to E lab rts 350 MeV. Since E lab < 2h 2 A 2 /M, M being the nucleon mass, 
and one requires Vi ow -k to reproduce the empirical phase shifts, a choice 
of A in the vicinity of 2 fm -1 seems to be appropriate. In our calculations 
we have used A = 2.1 fm -1 . 



3. Realistic Hartree-Fock Calculations 

As stated before, Vi ow -k is energy independent and well-behaved, so that 
one can use it directly to compute the bulk properties of 16 O in the frame- 
work of the self-consistent Hartree-Fock (HF) theory. The HF equations 
are solved using a basis of harmonic-oscillator wave- functions. The single- 
particle states \i) are expanded in a finite series of oscillator wave-functions 
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Figure 2. Binding energy per nucleon versus inverse r.m.s. mass radius for 16 O. The 
experimental data are from Rcfs. 15,16. 



\a): 



|i> = £C™» 



(6) 



The expansion coefficients C % a and the oscillator parameter fko are deter- 
mined solving self-consistently the HF equations 



£ 





(a\T\(3) +J2(<xJh\Vi ow - k \(3,j h ) 



a 



(7) 



Since we assume that 16 is a spherical nucleus, the \i)'s have good orbital 
and total angular momentum. The Coulomb correction is included exactly 
by taking the total two-body interaction as Vi ow -k plus the Coulomb po- 
tential. 

Once Eqs.(ll) have been solved, we use the HF wave-functions to calcu- 
late both the binding energy and r.m.s. mass radius. We also evaluate the 
second-order contribution to the above quantities within the framework of 
the Goldstone linked cluster perturbation expansion 15 . In Fig. 2 we report 
the results obtained using all the iViV potentials mentioned in the Intro- 
duction. It is worth noting that the differences between the various NN 
potentials are not very relevant, and that second-order corrections improve 
remarkably the agreement with experiment. The latter point confirms the 
need to go beyond a mean-field description of ground-state properties of 
finite nuclei 18 ' 19 - 20 ' 21 . 
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4. Summary and Conclusions 

In the present work, we have used a new technique to derive a low- 
momentum NN potential Vi ow -k from realistic potentials by integrating 
out the high-momentum modes. With this smooth Viow-k we have calcu- 
lated the bulk properties of 16 O in the framework of the Goldstone expan- 
sion using a self-consistent Hartree-Fock basis. The quality of the results, 
obtained using several modern realistic NN potentials, is quite good. We 
conclude that this new method, that it has also been successfully employed 
in shell- model calculations 12 ' 22 , may contribute to a better understanding 
of the role of realistic NN interactions in nuclear structure. 

Acknowledgments 

This work was supported in part by the U.S. DOE Grant No. DE-FG02- 
88ER40388, and the Italian Ministero dell'Istruzione, dell'Universita e della 
Ricerca (MIUR). 

References 

1. R. Machleidt, K. Holinde, and Ch. Elster, Phys. Rep. 149, 1 (1987). 

2. R. Machleidt, Phys. Rev. C 63, 024001 (2001). 

3. M. Lacombe, B. Loiseau, J. M. Richard, R. Vinh Mau, J. Cote, P. Pires, and 
R. de Tourreil, Phys. Rev. C 21, 861 (1980). 

4. V. G. J. Stoks, R. A. M. Klomp, C. P. F. Terheggen, and J. J. de Swart, Phys. 
Rev. C 49, 2950 (1994). 

5. D. R. Entem, R. Machleidt, and H. Witala, Phys. Rev. C 65, 064005 (2002). 

6. P. Lepage, nuc-th/9706029. 

7. D.B. Kaplan, M.J. Savage and M.B. Wise, Phys. Lett. B 424, 390 (1998). 

8. E. Epclbaoum, W. Glockle, and Ulf-G. Meissner, Nucl. Phys. A 637, 107 
(1998). 

9. E. Epclbaoum, W. Glockle, A. Kriiger, and Ulf-G. Meissner, Nucl. Phys. A 
645, 413 (1999). 

10. P. Bedaque et al. (eds.), Nuclear Physics with Effective Field Theory II 
(World Scientific, Singapore, 1999). 

11. S. Bogner, T.T.S. Kuo and L. Coraggio, Nucl. Phys. A 684, 432c (2001). 

12. Scott Bogner, T. T. S. Kuo, L. Coraggio, A. Covello, and N. Itaco, Phys. 
Rev. C 65, 051301(R) (2002). 

13. E. M. Krenciglowa and T. T. S. Kuo, Nucl. Phys. A bf 342, 454 (1980). 

14. F. Andreozzi, Phys. Rev. C 54, 684 (1996). 

15. J. Goldstone, Proc. Roy. Soc. (London) A 239, 267 (1957). 

16. G. Audi and A. H. Wapstra, Nucl. Phys. A 565, 1 (1993). 

17. A. Ozawa et al, Nucl. Phys. A 691, 599 (2001). 

18. K. Suzuki and R. Okamoto, Prog. Theor. Phys. 92, 1045 (1994). 



6, 2008 21:1 WSPC/Trim Size: 9in x 6in for Proceedings 



cortona2002 



7 

19. J. H. Heisenberg and B. Mihaila, Phys. Rev. C 59, 1440 (1999). 

20. H. Miither and A. Polls, Prog. Part. Nucl. Phys. 45, 243 (2000). 

21. A. Fabrocini, F. Arias de Saavedra, and G. Co, Phys. Rev. C 61, 0443002 
(2000). 

22. L. Coraggio, A. Covello, A. Gargano, N. Itaco, T. T. S. Kuo, D. R. Entem, 
and R. Machleidt, Phys. Rev. C 66, 021303(R) (2002). 



